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ABSTRACT 

Wc give an improved estimate of the detectability of gravitational waves from mag- 
netically confined mountains on accreting neutron stars. The improved estimate in- 
cludes the following effects for the first time: three-dimensional hydromagnetic ( "fast" ) 
relaxation, three-dimensional resistive ("slow") relaxation, realistic accreted masses 
^ 2 X 10~^Mq, (where the mountain is grown ab initio by injection), and ver- 
ification of the curvature rescaling transformation employed in previous work. Typ- 
ically, a mountain does not relax appreciably over the lifetime of a low-mass X-ray 
binary. The ellipticity reaches e « 2 x 10~^ for Ma ~ 2 x 10~^Mq. The gravita- 
tional wave spectrum for triaxial equilibria contains an additional line, which, al- 
though weak, provides valuable information about the mountain shape. We evalu- 
ate the detectability of magnetic mountains with Initial and Advanced LIGO. For 
a standard, coherent matched filter search, we find a signal-to-noise ratio of d = 
28{Ma/10-'^MQ){l + 5.5Ma/10-'*AfQ)-i(i:)/10kpc)-i(ro/14d)i/2 for Initial LIGO, 
where D is the distance and Tq is the observation time. From the nondetection of 
gravitational waves from low-mass X-ray binaries to date, and the wave strain limits 
implied by the spin frequency distribution of these objects (due to gravitational wave 
braking), we conclude that there are other, as yet unmodelled, physical effects that 
further reduce the quadrupole moment of a magnetic mountain, most notably sinking 
into the crust. 

Key words: accretion, accretion disks - stars: magnetic fields ~ stars: neutron - 
pulsars: general 



1 INTRODUCTION 

Accreting neutron stars in low-mass X-ray binaries 
(LMXBs) are promising sources of continuous gravita- 
tional waves (GWs). The signal from these emitters can 
be coherently integrated, so that the signal-to-noise ra- 
tio increases with the s quare root of the observation time 
ijjaranowski et al]|l998l ). Recent directed searches for GWs 
from the nearby X-ray source Sco X— 1 by the Laser Inter- 
ferometer Gravitational Wave Observatory (LIGO) set an 
upper bound on the gravitational wave strain of ho < 10~^^ 
l|Abbott et al.ll2007l ). 

LMXBs emit continuous gr avitational waves via a va- 
riety of physical mechanisms (jOwenl l2006l : lAbbott et al.l 
I2OO7I ): nonaxisymmetric elastic deformations of the neutron 
star crust, generated by temperature gradients (|Bildstenl 
1 19981 : lUshomirskv et al.l i2000l : iHaskell et all |2007| ) or inter- 
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nal toroidal magnetic fields (|Cutlei|[200^ ): r-modes, gen- 
erat ed by the Chand r asekhar-Friedmann-Sch u tz instabil- 
ity |Owen et all Il998l: lAndersson et al.1 Il999l : IStergioulaj 
l200a : iNavvar fc OwenI I2OO6I ): free precession, excited 
by i n ternal or accret i on to r ques (jjones fc AnderssonI 
20021: IVan Den Broec^ l2005l : iPavne fc MelatosI bOOGaT 



Chung et al.1 |2008|) 



tains |Pavne fc Melatoi 



and mag n etically 
20041 



confined moun- 
Melatos fc Pavnel l2005l : 



iPavne fc MelatosI l2006al : I Vigelius fc Melatodl2008al 



In the latter mechanism, accreting plasma accumu- 
lates at the magnetic poles and spreads equatorwards. The 
frozen-in magnetic field is carried along with the spread- 
ing plasma and is therefore compressed, to the point where 
magnetic tension counterbalances the latitudinal pressure 
gradient. This e quilibrium configuration is termed a mag- 
netic mountain (|Pavne fc Melato3 |2004| ) . During the pro- 
cess, the magnetic dipole moment of the star decreas es with 
accre t ed mass, consistent with o b servational data (IZhand 
1 19981 : iMelatos fc PhinnevI I2OOII : IPavne fc MelatosI 12004 : 
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IZhang fc Koiimal 120061 '). The distorted magnetic field can 
also act as a thermal barrier between th e hemispheres, af- 
fecting the physics of type I X-ray bursts (jPavne fc Melatoa 
l2006bl i. 

In this paper, we draw together the latest analytic 
and numerical modeUing of magnetic mountains in LMXBs 
jVigehus fc Melatod l2008al lb[) to compute rigorously the 
signal-to-noise ratio of these sources for a coherent search 
with LIGO. To this end, we make extensive use of previously 
published results. Hydromagnetic equilibria of magnetically 
confine d mountains were computed analytically and numeri- 
cally bv lPavne fc Melato s (,2Q04) ■ While t hese configurations 
are st able to axisymmetric perturbations (|Pavne fc Melato^ 
I2OO7I ). the (nonaxisymmetric) undulating submode of the 
three-dimensional Parker instability induces a recon figura- 
tion of the magnetic field (|Vigelius fc Melatos|[2008al ). How- 
ever, the line-tying boundary condition at the stellar sur- 
face prevents the mountain from being disrupted and the 
saturation state of the instability still exhibits a substan- 
tial quadrupole moment wit h a high degree of axisymmetry. 
IVigehus fc MelatosI (|2008bh extend the analysis to include 
resistive effects and demonstrate that the mountain decays 
on the diffusion time scale. No evidence for resistive insta- 
bilities that grow on a short time scale is found. 

In their stability a nalysis , JPavne fc MelatosI (|2006al ) 
and I Vigelius fc MelatosI (|2008al lb[) numerically computed an 
equilibrium configuration with a particular value for the ac- 
creted mass. Ma, subsequently loaded this equilibrium into 
a magnetohydrodynamic (MHD) solver and evolved it. In 
this article, we solve the initial-value MHD problem for 
the first time by injecting plasma into an initially dipo- 
lar field. This method allows us t o independently validate 
and extend the scalings giv en bv IPavne fc MelatosI (|2004l ) 
and iMelatos fc Pavnd (|2005l 'l , which were calculated analyt- 
ically in the small- Ma approximation. Taking into account 
resistive effects and three-dimensional reconfiguration of the 
mountain, we present easily applicable formulas to compute 
the mass ellipticity for a given Ma and give improved esti- 
mates on the strength of the gravitational wave emission. 
In particular, we investigate how the (small) degree of non- 
axisymmetry changes the gravitational wave spectrum and 
how gravitational-wave spectrometry can be used to obtain 
valuable information about the underlying field configura - 
tion. This analysis ties in with jPavne fc Melatoa (|2006al l 
who consider axisymmetric mountains. 

The paper is organised as follows. We describe quan- 
titatively the physics of the mass quadrupole moment of a 
magnetic mountain in section[2l including three-dimensional 
force balance, slow (resistive) and fast (hydromagnetic) re- 
laxation, realistic accreted masses, and the influence of stel- 
lar curvature. Taking into account these effects, we give a 
recipe to compute the quadrupole moment as a function 
of accreted Ma in section [3] and hence estimate the strength 
and detectability of the GW signal. The frequency spectrum 
of the signal is calculated in section (4] We discuss our results 
in the context of past and future LIGO searches in section 



2 MASS QUADRUPOLE MOMENT OF A 
MAGNETIC MOUNTAIN 

In the context of gravitational radiation, the key property 
of a magnetic mountain is its mass quadrupole moment, 

Qij = J d^x {Sx'ix'j - r''^Sij)p(x.'), (1) 

where p denotes the plasma density. We aim to calculate Qij 
as a function of the accreted mass Ada ■ It is useful to mea- 
sure Ma in units of the critical mass Mc = GMtB^Rl /{8ci), 
where M* and _R* are the stellar mass and radius, £?» = 10^^ 
G is the initial magnetic field, and Cs = 10* cm s~^ is 
the isothermal sound speed. For Ma > Mc, the magnetic 
dipole moment p decrea ses with Ma', for Ma < Mc, fJ, is 
approximately constant (jPavne fc Melato3l2004i ). In a typ- 
ical LMXB, we have 10"'' < Mc/Mq < 10"* and 10"^ < 
Ma/Me < 10" ^ that is, magnetic burial distorts the field 
dramatically. 

The three-di m ensiona l equilibria computed by 
IVigelius fc MelatosI (|2008al ) deviate from axisymmetry 
(with respect to the magnetic axis) by less than 0.1 per 
cent in the mass quadrupole moment. Hence, we frequently 
employ the axisymmetric mass ellipticity, defined as 

6=-^ / dfdrpr* sin ^(Scos^e - 1), (2) 

to describe the mountain as a biaxial ellipsoid, where J^j = 
2M,7?,/5 denotes the unperturbed moment of inertia. 

In order to give a reliable estimate of the mag- 
nitude of e and hence the strength of the gravitata- 
tional wave sig nal, we take into account the effects of 
hydromagnetic jVigelius fc MelatosI l2008al ) and resistive 
(| Vigelius fc MelatosI l2008bh relaxation and consider the 
high-Ma limit as well as the influence of the neutron star 
curvature. For the convenience of the reader, we repeat pre- 
viously published results in the first three subsections. We 
compute Qij for the three-dimensional equilibrium state in 
section [2711 We then examine how the mountain responds to 
hydromagnetic and resistive relaxation in sections 12.21 and 
12.31 respectively. Section 12.41 explains how to build moun- 
tains with realistic values of A'la , and the effect of curvature 
downscaling is quantified in section 12.51 The results in the 
last two subsections are new and have not been published 
elsewhere. 



2.1 Three-dimensional equilibrium 

In a magnetic mountain at equilibrium, the pressure gra- 
dient balances the gravitational and Lorentz forces. Start- 
ing with a centred magnetic dipole before accretion be- 
gms, IPavne fc MelatosI (|2004l ') computed the unique, self- 
consistent, Grad-Shafranov equilibrium that satisfies force 
balance, while simultaneously respecting the flux-freezing 
constraint of ideal magnetohydrodynamics (MHD). 

The top-left panel of Fig. [1] displays an axisymmetric 
equilibrium for Ma = Mc. The mountain (orange isosur- 
face) is confined to the magnetic pole by the tension of the 
distorted magnetic field (blue and green curves). Blue and 
green field lines are drawn starting from the pole and equa- 
tor, respectively. The region where the magnetic pressure is 
greatest (at 8 « 7r/4 in Fig. (1^) is termed the equatorial 
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Figure 1. Density isosurfacc (orange) and magnetic field lines (blue and green) for a magnetic mountain, (a) Axisymmetric, idcal-MHD 
configuration with Ma = 1-2 X 10~*Mq. (b) Nonaxisymmetric, ideal-MHD configuration with Ma = 1.2 X 1Q~'^Mq. (c) Snapshot of the 
resistive evolution at t = IOta, with Ma = 1-2 X IO^^Mq. (d) A mountain grown by injection with Ma = 1.9 x 10~^Mq. The mountain 
is defined by the orange isosurface p{r, 9, if)) = O.Spmax, where pmax = 2.0 X 10® g cm~'^ is reached at x = (r — Rt)/ho = 0.9 X 10~^ and 
6 = 0.01 in the axisymmetric model. In order to aid the reader, the altitude scales are magnified five-fold in all panels. The foot points 
of the blue field lines touch the stellar surface, while green field lines are traced starting from the equator. 



magnetic belt. Here, B is ~ 16 times higher than at the 
pole. 



However, an axisymmetric analysis neglects important 
toroidal modes. When we load the axisymmetric equilib- 
rium in Fig. [1^ into the ideal-MHD code zeus-mp, we ob- 
serve that it is unstable to the undulating submode of the 
three-dimensional Parker instability, which reconfigures the 
hydr omagnetic structure by gro wing the toroidal magnetic 
field (|Vigelius fc MelatosI [2008al ') . Thanks to line tying, the 
instability is not disruptive; its saturation state (Fig.[T]3) still 
confines matter at the magnetic pole, and the mountain is 
still present even though e is reduced to ~ 60 per cent of its 
original value. The evolution of Qij for this model (Fig. \^) 
is displayed in Fig. [S] The diagonal components, describing 
the axisymmetric distortions, decrease by « 40 per cent, on 
a time-scale of « 2ta- Here, ta = (7ri?,p^'''^/_B)niin denotes 
the pole-equator crossing time for transverse Alfven waves; 
nRtp^^^ / B is smallest close to the stellar surface. For a re- 
alistic star with Ma = Mc, the Alfven time-scale evaluates 
to TA = 5.1 X 10"^ s. 



2.2 Fast, hydromagnetic relaxation 

A magnetic mountain performs global hydromagnetic os- 
cillations when perturbed, but it remains intact. This un- 
expected outcome can be ascribed to two factors: (i) the 
mountain is already the saturation state of the nonlinear 
Parker instability, and (ii) the line-tying at the stellar sur- 
face su ppresses important localised modes, e.g. interchange 
modes. IPavne fc Melatoi (|2007l ) found all mountains with 
Ma ^ 6 X 1O~*M0 to be marginally stable. 

The mountain quickly tends to an almost axisymmet- 
ric state {\Qi2/Q:i3\ < in Fig. This high degree 
of axisymmetry considerably simplifies the computation of 
the amplitude of the gravitational wave strain (see section 
I3.1|l . (Note that Qij is defined relative to the magnetic 
axis, which is inclined with respect to the rotation axis). 
The reconfiguration is accompanied by global, nonaxisym- 
metric, MHD oscillations (top panel of Fig. [SJ. Although 
they die away in this numerical experiment, global oscilla- 
tions can be continuously excited in reality (e.g. by accre- 
tion torques), modifying the gravitati onal wave spectrum 
(see section l4)). |Pavne fc MelatosI (|2007h identified two dom- 
inant modes: a short-period sound mode, with a frequency 
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considerable uncertainties about the exact value of cr, such 





t/r. 

Figure 2. Quadrupolc moments for the nonaxisymmetric config- 
uration (Fig. [T})), normalised to Q33 = 1.30 X 10^^ g cm^ of the 
axisymmetric model, as a function of time, in units of the Alfven 
time, ta = 5.1 x 10"'^ s. The system develops an appreciable 
nonaxisymmetry, characterised by the off-diagonal element Qi2> 
during the relaxation phase before settling into a nearly axisym- 
metric state. 



/s/kHz = 1.4x 10'^(cs/10*cms"^) (independent of Ma), and 
a longer period Alfven oscillation, which can be fitted by 
/a = 17(Ma/M,) Hz. 

(ILitwin et al ] |200ll ) demonstrated that an accretion col- 
umn that is magnetically confined to the polar cap is sus- 
ceptible to the ideal-MHD ballooning instability. However, 
our simulations do not exhibit any evidence for a growing 
instability. Physically, this is because the compressed equa- 
torial magnetic field stabilises the lateral motions involved 
in such an instability. Indeed, in Fig. [S] we demonstrate 
that a magnetically confined mountain is only susceptible 
to a growing ballooning mode when the back-reaction of 
the magnetic belt is neglected. We perform an axisymmet- 
ric simulation with Ma = Mc and ^ ^ ^ tt/S, where the 
outer ^-boundary (at 9 = it/8) is set to outflow. The moun- 
tain is clearly disrupted on the Alfven timescale by the bal- 
looning mode, as in Fig. [S] the magnetic field and frozen-in 
plasma slide sideways through the outfiow bounda ry. This is 
consistent with the findings of iLitwin et al.l (|200ll ). who im- 
posed boundary conditions equivalent to outfiow, therefore 
neglecting the stabilizing efi^ect of the equatorial magnetic 
belt. Furthermore, our growing simulations do not show any 
evidence for an instability during the early stages of accre- 
tion (when Ma < Mc). 



2.3 Slow, resistive relaxation 

A magnetic mountain relaxes resistively over a long time- 
scale, which is set by the ohmic diffusion rate across the 
steepest magnetic gradients in the mountain. Resistive re- 
laxation reduces Qij. Simple estimates suggest that the 
effect can be neglected as long as Ma < 1O~"'M0, as- 
suming a homogeneous conductivity, a — 7.7 x 10^'' s^^ 
iMelatos fc Pavnell2005l ). Here, we assume that a is dom- 
inated by electron-phonon scattering, with a crustal tem- 
perature T = 10^ K an d a characteristic pla sma density 
of /9 = 5 X 10" g cm"^ (|Cumming et al.l 120041 '). There are 



as the value of the impurity parameter JSchatz et al.| 19991: 
Gumming et~al] I2001I . |2004| : Ijonesi |2004| : iPons fc Geppert 



2007 ) and the composition o f the crust (jGumming et al 
20041 : IChamel fc Haenseill2008l 'l. Furthermore, a is in reality 



a function of the lo c ation t hrough T and p. In keeping with 
IVigeUus fc MelatosI (|2008bh . we treat the electrical conduc- 
tivity as a fiducial parameter and note that the time-scale of 
resistive relaxation scales with a. An inhomogeneous con- 
ductivity will be considered in a forthcoming paper. 

Transient resistive instabilities , like a global te aring 
mode or local gravitational mode (|Furth et al.l 1 19631 ) , are 
known to evolve on time-scales as short as ^ (tota)^^^ ~ 
0.3 yr, where td denotes the characteristic diffusion time- 
scale. Note that td implicitly depends on the position 
through B and p and has a minimum close to the sur- 
face. The definition of the characteristic diffusion timescale 
is hence somewhat arbitrary. As justified in section [2m we 
pick the minimum value td ~ (47r(ji3/c^| V^B|^ . . Insta- 
bilities grow in magnetic neutral sheets [creat ed by the un- 
dulating submode of the Parker instability; see lHanasz et al.l 
((2002)] or regions of hig h magnetic shear. 



Vigelius fc MelatosI (|2008bl ) tested numerically whether 



neutral sheets can grow to disrupt the mountain on short 
time-scales. They evolved the three-dimensional equilibrium 
in Fig. [TJj, and similar states for other values of Ma, in 
ZEUS-MP, extended to treat ohmic diffusion. The results are 
reported in Fig. [l] which depicts the evolution of e for differ- 
ent Lundquist numbers Lu — td/ta. For each value of Lu, 
e{t) is plotted as a function of time, measured in units of the 
Alfven tim^B ta = 2.5 x 10~^ s (left panel) and the diffusion 
time Td (right panel), which differs for each model. The mod- 
els with a realistic resistivity {Lu = lO^**, solid curve) and 
Lu = 10'' (dotted) do not exhibit any change in e over the 
simulation time. For Lu = 10^ (dashed), we note a decrease 
of 21 per cent over O.Itd. For Lu = lO"'^ (dash-dotted), 
the mountain relaxes substantially; e drops by 90 per cent 
over the diffusion time-scale. For Lu = 10~^ (dash-triple- 
dotted), the mountain immediately slips through the mag- 
netic field lines and falls freely towards the equator, where 
it is reflected by the boundary surface, causing e to oscillate. 
The latter case, in particular, is of academic interest only, as 
far as its application in LMXBs is concerned. Magnetic neu- 
tral sheets are found in the toroidal plane, where the plasma 
density and magnetic field strength are low. Reconnection 
occurs locally in these regions, smoothing toroidal gradients. 
Ohmic dissipation therefore tends to restore axisymmetry. 

Importantly, the mountain relaxes globally on the dif- 
fusion time-scale, td, which greatly exceeds the accretion 
time-scale race, with e falling to e~^ of its initial value af- 
ter 34rD. In practice, this means that the three-dimensional 
saturation state of the Parker instability does not relax re- 
sistively until ~ 10^ yr elapse, at least for the examined 
moun tains with Ma S, W^^Mq . In the analytic small- Ma 
limit, iMelatos fc Pavnd (|2005h found the minimum accreted 
mass for ohmic diffusion to arrest mountain growth to be 
Md ^ 10~^Mq, provided the star accretes at the Eddington 



^ In Fig.|4] TA is the characteristic Alfven time-scale for the three- 
dimensional equilibrium. It is half the characteristic Alfven time- 
scale of the axisymmetric configuration used in section 12.11 



Detectability of gravitational waves from magnetic mountains 5 




Figure 3. Meridional section of a mountain with Ma = Mc and outflow boundary conditions at 8 = n/8 at t = (left panel) and 
t = 0.9ta (right panel). The mountain is susceptible to an ideal-MHD ballooning instability and disrupts over the Alfven timcscale. 
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Figure 4. Evolution of mass ellipticity e for different values of the Lundquist number Lu = to/ta = 10^^, 10"^, 10"^, 10~^, 10"'^ (solid, 
dotted, dashed, dash-dotted, dash-triple-dotted), curves from top to bottom, in units of the Alfven time ta = 2.5 X 10~^^ s (left panel) 
and the respective diffusion times td (right panel). The solid and dotted curves nearly overlap. The magnetic mountain relaxes at the 
time-scale tq, with e falling to e^^ of its initial value after 34r]3. Two of the curves are hard to see in the right panel: Lu = 10^* (solid) 
stops at f = td and Lu = 10^ (dotted) stops at t = lOro. 



rate (such that race = M^/Ma ~ 10^ yr) and has a crustal 
temperature of T = 10* K. Our numerical result, indepen- 
dently validated in section 12.41 raises Md by three orders of 
magnitude. 

In addition to resis t ive eff ects in the accreted plasma, 
iKonar fc Bhattacharval (| 19971 ) explored the influence of 
ohmic dissipation on the crustal magnetic field. They iden- 
tified three competing mechanisms: (i) the accretion fiow 
advects current into lower layers with higher p and conse- 
quently higher a; (ii) the crust is heated by advection, de- 
creasing a; and (iii) the current is squeezed into the inner 



layers, decreasing td- As a result, the crustal field decays 
rapidly before freezing at a residual surface value, whose 
magnitude increases with Mg [by reducing t he duration of 
the rapid decay phase, see also [Roman il(|l990l)]. We defer the 
study of sinking and stratified a [compare, e.g., (|jahan-Miril 
I2OO0I : IChoudhuri fc KonaillioO^ )] to future work. 
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2.4 Realistic accreted mass 

Numerical obstacles, like steep gradients, and physical ob- 
stacles, like magnetic bubble formation, interfere with the 
task of modelling magnetic mountains for realistic val- 
ues of Mg. The i terati ve numerical scheme employed by 
IPavne fc MelatosI (|2004 ) to compute Grad-Shafranov equi- 
libria converges poorly for Ma > lO~*Af0. A bootstrap- 
ping algorithm can be used to quasistatically fatten a G rad- 
Shafranov equilibrium ten-fold (|Pavne fc Melato3l2007l 'l , but 
the results have not yet been verified against a reliable 
numerical solution of the full initial-value problem (boot- 
strapping converges quite violently), and anyway, bootstrap- 
ping works up to Ma ~ 10~^Mq at most. In short, a self- 
consistent configuration w ith a realistic amoun t of accreted 
matter, e.g. Ma ~ O.IM© l|Burderi et al.lll999l ). is yet to be 
achieved. 

We try to overcome this restriction in this subsection 
by growing a magnetic mountain from scratch, by inject- 
ing plasma at the inner boundary r = R, into an initially 
dipolar background field. The injection speed is chosen to be 
less than the gravitational escape speed from the simulation 
box, while the density is chosen to give ta <C Tacc, such that 
the system passes through a sequence of quasistatic equilib- 
ria, yet Tacc is short enough to keep the simulation runtime 
reasonable for Ma ~ 10~^ Mp). This approach diffe rs from 
the bootstrapping algorithm (|Pavne fc Melato^l2007l ) in two 
ways: (i) instead of relying on a Grad-Shafranov equilibrium 
as the starting point, we solve the full initial- value problem; 
and (ii) we inject plasma at the r = R, boundary, thereby 
circumventing the artificial field line pinning at the outer 
boundary that stems from the inflow boundary condition. 
Implementation details are provided in Appendix 1X1 

An axisymmetric grown mountain with Ma = 1.9 x 
10~^Mq and & = 3 (hemispheric-polar magnetic flux ra- 
tio; see Appendix IX)) is displayed in panel (d) of Fig.[T] The 
mountain isosurface covers the whole star. The base density, 
at {x,6) = (10"^, 0.012) is fifty times higher than for Ma = 
1.2 X IO^'^Mq [panel (a)]. At first glance, the magnetic field 
configuration looks entirely different: instead of pointing ra- 
dially outward, all field lines are closed loops. This ostensible 
difference is due to the boundary condition dH/dr = at 
r = Rm, enforced i n the growing simu l ations (see Appendix 

cL Be = in IPavne fc MelatosI (|2004 l. However, ap- 
pearances are a bit misleading: the all-important equatorial 
belt, where the magnetic field is highly distorted and most 
intense, is clearly visible in Fig.[T]i, just as much as Figs. [T^ 
-d]:. Near the pole, at {x,e) = (10"^, 0.012), B in model 
d is comparable to B in model a. At the magnetic equa- 
tor, {x, 9) — (10"'^, 1.4), _B is ~ 25 times higher in model d 
than in model a. Naturally, the magnetic tension required to 
counterbalance the hydrostatic pressure is greater in model 
d. 

Fig. [5] shows the evolution of e as a function of Ma up 
to Ma ^ 1.9 X lO"^Af0. Plotted are models with zero resis- 
tivity (solid curv e ), a r ealistic resistivity r] — 1.3 x 10~^^ s 
[Gumming et al.l (|2004l ') . plus symbols], and an artificially 
high resistivity, 77 = 9.2 x 10^^^ s (i.e. Lu — td/ta = 
10^ dashed curve). We also fit a functional depen dence 
similar to that propos ed by IShibazaki et al. 1 (|1989D and 
iMelatos fc Pavng (|2005l ) to the plus symbols, finding 
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Figure 5. Mass cUipticity e as a function of Ada- The solid 
curves show a model with zero resistivity (solid curve), realis- 
tic resistivity rj = 1.3 X 10^^^ s (plus symbols), and an artificially 
high resistivity rj = 9.2 X 10^^^ s (dashed curve). The dotted 
curve represents a fit to the (slightly modified) formula given in 
iMelatos fc Pavnd 1I2OO5I) : see equation l|3). 

« -I -I Ma f, , Ma V 

10^ = ^^irao (,' + '-'iraej ■ 

The fitting formula is plotted as a dotted curve in Fig. [S] 

Fig.[5]and equation ([Sj are indispensable tools for calcu- 
lating the gravitational wave emission, summarising as they 
do our most up-to-date modelling of resistive relaxation and 
large (i.e. realistic) values of Ma- For Ma > 0.5 x 10"'' Mq, 
e is expected and found to follow the shape of the small- 
Ma analytic solution, valid for Ma <g: Mc = 1.2 x 10"* [see 
Appendix in iPavne fc Melatoi (|2004 )]. The deviation visi- 
ble at Ma = lO"^Af0 can be attributed to the fact that 
the configuration has not yet had time to equilibrate. In 
this regime, material piles up in a polar flux tube of ap- 
proximately constant cross-sectional area, yielding e oc Ma- 
For Ma > 10~^A-/q, the hydrostatic pressure overcomes 
the magnetic tension and the mountain spreads towards 
the equator, triggering global MHD oscill ations. These lat- 
eral o scillations, observed numerically bv lPavne fc MelatosI 
(|2007h . compress the magnetic fleld once per cycle, increas- 
ing |V^B|, promoting lateral diffusion of the mountain, and 
reducing e. For high 77 (dashed curve), the oscillation in e is 
clearly visible in Fig. [Sj It is also visible, albeit less promi- 
nently, for realistic rj (plus symbols in Fig. [5] Hence the 
fairest way to interpret the e-Ma trend in Fig. [5] is that e 
saturates for Ma ^ Mc (fiat underlying trend) with an os- 
cillation superposed. 

In a realistic accretion scenario, given Ma, we expect 
to find one particular value of Ma (and hence e) at which 
the system attains a steady state, where the mass diffuses 
through fiux surfaces at a rate which is exactly replenished 
by Ma ■ A full parameter study to examine the relation e{Ma) 
is outside the scope of this a rticle. We simply note, following 
IVigelius fc MelatosI (|2008bh . that magnetic mountains are 
resistively stable over Tacc or A/a ^5 1.2xlO^''Af0 (neglecting 
mass inflow). 

Can short-timescale instabilities that are absent in the 
equilibrium configuration grow during the early stages of ac- 
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cretion, i.e. for small Ma? In this case, high-Ma equilibria 
could never be reached and magnetic mountains would not 
emit detectable gravitational radiation. However, our grow- 
ing simulations (Fig. [Sjl show no evidence for ideal-MHD or 
resistive instabilities during the low- A/a stage of accretion. 



2.5 Realistic curvature 

In general, the characteristic length-scale for radial gradi- 
ents (ho) is much smaller than the length-scale for latitu- 
dinal gradients (R*), creating numerical difficulties. How- 
ever, in the small- Ma l i mit, it can be shown analytically 
jPavne fc Melatod |2004 bOOeal ) that the structure of the 
magnetic mountain depends on R, and A'l, through the 
combination ho oc R^ /M* , not separately. We therefore arti- 
ficially reduce R, and M*, while keeping ho fixed, to render 
the problem tractable computationally. It is vital to bear in 
mind that invariance of the equilibrium structure under this 
curvature rescaling does not imply invariance of the dynam- 
ical behaviour, nor is the scaling necessarily applicable at 
large Ma- 

A standard neutron star has M* = 1.4M0, i?* = 10^ 
cm, Bt — 10^^ G, and Cs — 10* cm s~^. We rescale the star 
to M» = 1.0 X IO'^Mq and K = 2.7 x 10^ cm, reducing 
a = R*/ho to 50 while keeping it large. We then calcu- 
late the mountain structure, and hence e numerically. We 
upscale Ma back to a realistic star, using the scaling rela- 
tion for Mc, the critical accreted mass above which the star's 
m agnetic moment starts to change, defined by equation (30) 
of lPavne fc MelatosI (|2004 l: 



Me 
M« 



(4) 



Furthermore, we use the analytic result e oc 
l|Melatos fc Pavne|[2005l '). valid for Ma < Mc, to upscale e. 
In order to verify the fairness of this procedure, we perform 
runs for a = 100 and a — 500. If the curvature rescaling 
is fair, these runs should obey e cx a^. Computational costs 
limit us to achieving Ma = 0.34 x 10"'* Mq. We compute e, 
upscaled to a realistic star, for the maximum Ma for both 
runs, finding a relative deviation between the simulation out- 
put and the predicted scaling of < 10"^ per cent. 



3 GRAVITATIONAL WAVE STRAIN 
3.1 Ellipticity 

The neutron star and the piled up matter at the mag- 
netic pole can be modelled approximately as a rigid, bi- 
axial top, which is symmetric about the pre-accretion mag- 
netic axis. Biaxial equilibria are of course unstable to the 
toroidal Parker mode and reconfigure i nto a triaxial equi- 
librium , as described in section 1^?^ and [Vigelius fc MelatosI 
(|2008al ). However, the ultimate deviations from axisymme- 
try are small, less than 0.1 per cent, according to Fig. |2l 
and can be neglected in a first analysis. In this case, we 
can rewrite equation ([2]) as e — j/a — 7i|//i, where Ii and 
/a denote the moments of inertia with respect to the star's 
principal axes. It is important to keep in mind that we are 
dealing with a prolate spheroid, with Ji > J3; the conse- 
quences for the long-term rotational evolution are explored 



in section |4~2] Frequently, authors omit taking the absolute 
value in the definition of e resulting in a negative ellipticity 
for our case. 

In general, the mountain axis is tilted with respect 
to the rotation axis. Thus the neutron star precesses 
freely, generating gravitational waves at /, and 2/,, where 
/* = J is the star's spin frequency. For a biaxial 
star, the wave strains in orthogonal (-1- and x) polarisa- 
tions can be written as (jZimmermann fc Szedenit^ Il979l : 
IJaranowski et al"]| 19981 ') 



h+{t) = -/iosin2Ssin2icos<I>(f) 
8 

-l-i/io sin^ 9{l + cos^ i) cos 2'l>(t), 

hx{t) = i/to sin 20 sin i sin $(t) 
+ho sin^ 6 cos i sin 2$ (t ) . 



(5) 



(6) 



Here, 6 denotes the wobble angle (between the total angular 
momentum J and principal axis of inertia es), i is the incli- 
nation angle (between J and the line of sight, drawn from 
the star to the solar system barycenter), ho is a character- 
istic amplitude. 



ho 



c4 D 



(7) 



(where I is the moment of inertia and D is the distance to 
the source), and <I>(t) is the phase, including Doppler terms; 
for a source at rest relative to the observer, we can write 
$(f) = 27r/.t-H$o. 

Upon combining all the results in section |3] principally 
equation Q, Fig. [S] and the multiplication ~ 0.6 for con- 
verting e from two- to three-dimensional equilibria (section 
I2.1|l . we arrive at the following approximate formula for the 
ellipticity in the absence of resistivity: 



10- 



6.82 



Ma 



10-*M, 



© 



1 + 5.5 



Ma 



10-*Me 



(8) 



Resistive relaxation arrests the growth of the mountain 
(and hence e) at a value of Ma, denoted Md, which depends 
on Afa (see section I2.4|l . Conservatively, we conclude from 
the results in section that a realistic resistivity does not 
relax a mountain with Ma = 1.2 x 10~*Mq over the accre- 
tion time-scale. This implies e ^5 3.6 x 10^**, but the true 
value of e is expected to be much lower than the upper 
bound. 



3.2 Signal-to- noise ratio 

The signal x{t) read out at the detector port is buried in 
noise. We assume here that the noise n{t) is additive, sta- 
tionary, and Gaussian, with {n{t)} — 0. Then we can write 
IJaranowski et al]| 19981 ) 



x{t) = h{t) + n(t), 
where 

h{t) = F+{t)h+{t) + Fx{t)hAt) 



(9) 



(10) 



contains the beam pattern functions F+ and Fx , which en- 
code the diurnal motion of the Earth. 
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By averaging over sky position, inclination, and polari- 
sation, the signal-to-noise ratio can be expressed as a func- 
tion of the wobble angle, ho, and the observation time To, 
viz. 



(4) 



and 



hlTo sin^ 261 
1005h(/,) ' 



sin* e 



(11) 



(12) 



25S42/.) ' 

assuming that the interferometer arms are perpendicular. 
Here, Sh{f) denotes the one-sided spectral noise density of 
the detector. 

It is sometimes desirable to average over wobble angle 
as well, i n the absence of knowledg e about a specific object. 
Following iPavne fc Melatod (|2006al 'l. we average m a manner 
that is unbiased towards small 9, viz. {) ~ d(cosS); cf. 
iThornd (|l987l 'l. The final resuh is 



2hlTo 
375 



16 



(13) 
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3.3 LIGO detectability 

In Fig. [S] we plot d as a function of accreted mass and 
wobble angle 6 for a standard pulsar with / = 500 Hz, lo- 
cated at a distance D = 10 kpc, assuming a coherent in- 
tegration time of To = 14 d. The inputs are equation ((T)), 
equation ((8]), equation (|ll|l . and equation l|12p . We neglect 
resistive relaxation in order to obtain an upper limit on d; 
the results of section 1^751 and section [5TT] suggest Ohmic dif- 
fusion saturates d for Ma S5 1O~*M0. The solid and dashed 
curves in Fig. [6] refer to the current and planned sensitiv- 
ities of Initial and Advanced LIGO respectively, published 
in LIGO science requirement document and the Advanced 
LIGO proposaffj. We find that a perfectly radiating moun- 
tain with Ala ^ 1O~*M0 is barely detectable with Initial 
LIGO (d ~ 1) and firmly detectable with Advanced LIGO 
(rf ~ 10). As Ma increases, d saturates at ~ 10 for LIGO and 
~ 100 Advanced LIGO. One must remember, though, that 
resistive relaxation reduces these values for Ma > 1O~*M0. 

An alternative way to estimate the detectability of 
the signal is to compare the characteristic wave strain 
ho versus the statistical threshold hth- A signal is de- 
tected with a false alarm rate of 1 per cent and a false 
dismissal rate of 10 per cent when hp e x ceeds hth ~ 
11.4[ gh(2/.>)]^/^rn~^ iJaranowski et al.l Il998l : lAbbott et al] 
I2OO4I '). Fig. [7] displays ho as a function of Ma for the same 
pulsar as in Fig. [5] The diagonal lines give ho for Ma/M^ = 
10"^ 10"*, 10"^ 10"^ 10"^ lO-*, lO-^, for & Waxlal (solid 
lines) and a triaxial (dashed lines) mountain. Also plotted 
as solid curves are /ith for Initial LIGO and Advanced LIGO 
for To = 14 d, averaged over 9 and i. We exclude the re- 
gion Ma > 1O~*M0 (light shaded region), in which resis- 
tive relaxation prevents further growth of e, and the region 
2/ > 1.4 kHz (dark shaded region), because no accreting 
millisecond pulsars have been discovered yet in this band 
(they may be in the future). 

Fig. [7] suggests that there is a small region 400 < 
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Figure 7. Amplitude ho of the gravitational wave signal for 
Ma/MQ = 10"^, 10-*, 10"^, 10"^, lO"'^, 10-"*, 10"^, for axisym- 
metric (solid diagonal lines, bottom to top) and nonaxisymmet- 
ric (dashed lines) mountains. The sensitivities of Initial and Ad- 
vanced LIGO, assuming 14 days coherent integration, are also 
plotted (upper and lower curves respectively). The growth of the 
mountain is arrested for Ma 1O~*M0 (light shaded region) by 
resistive relaxation llMelatos fc Pavnd ISooil : IVigel ms&Melatoa 
l2008alV The right-hand edge is excluded at present because no 
accreting millisecond pulsars have been discovered with 2/* > 1.4 
kHz (dark shaded region). 



/,/Hz < 1400 for which magnetic mountains are, in prin- 
ciple, detectable with initial LIGO. However, there are sev- 
eral physical mechanisms not yet considered in our mod- 
elling, most notably the sinking of the mountain into the 
crust, which act to reduce the gravitational wave signal. We 
discuss these mechanics further in section [5] 



4 GRAVITATIONAL WAVE SPECTRUM 

In this section, we investigate the gravitational wave spec- 
trum in more detail. Two effects modify the spectrum away 
from its simplest form (delta functions at /* and 2/*): global 
MHD oscillations, and triaxiality. We study the former with- 
out the latter in section HTT1 and vice versa in section to 
isolate the physics of the two effects. 



4.1 Biaxial, vibrating mountain 

A magnetically confined mountain oscillates when plucked, 
e.g. by starquakes or fluctuations in the accretion torque 
(|Pavne fc Melatodl2006al : I Vigelius fc Melatosll2008ah . These 
global hydromagnetic modes appear in the gravitational 
wave spectrum as acoustic and Alfvenic sidebands beside 
the two main peaks at /* and 2/». 

The top panels of Fig. |8] display the truncated Fourier 
transform. 
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Figure 6. Signal-to-noise ratio d for 10~^ ^ Ma/MQ ^ 10~^ for the nonaxisymmetric equilibrium for Initial LIGO (solid curve) and 
Advanced LIGO (dashed curve) after 14 days of coherent integration. The wobble angles are: 6 = it/A (top left), 9 = tt/S (top right), 
6 = it/2 (bottom left), unbiased average over ^ 6 ^ 7r/2 (bottom right). 



Hf) = jrj\te''^f'h{t), (14) 

of the wave strains in the + and x polarisations, given 
by equations ((Sjl and ((6|, for a standard star with Ma = 
1.2 X IQ-^Mq at D = 10 kpc, with 6 = n/3 and i = n/3, 
assuming an integration time of To — 14 days. For com- 
parison, we also plot hth- It is important to bear in mind 
that, for detection, one uses the combined signal power in 
both polarisations, even when the peaks in each single po- 
larisation remain under the threshold. The wobble angle is 
chosen artificially large to illustrate the effect of precession. 
For 6 = 7r/3, the two peaks at /, and 2/, have similar 
strengths. The width of the peak, the lopsidedness, and the 
signal power in between are numerical artifacts caused by 
the discrete Fourier transform; they remain even when we 
input the unmodulated e into equation (1 7)). For 9 = 0.1, a 
more realistic choice [cf. PSR B1828-11. iLin^ (|2003l )]. the 
peak at 2/* only reaches 28 (11) per cent of the /, peak in 
the + ( X ) polarisation (bottom panels of Fig. [Sjl . When the 
total signal power is distributed into two peaks of similar 
height, coincidence experiments at two frequencies become 
possible. 

The acoustic and Alfvenic sidebands produced by the 



global oscillations are very hard to see in Fig.[8l due to strong 
numerical damping. In a realistic neutron star, the Alfven 
mode may be perpetually re-excited (e.g. by fluctuations 
in the accretion torque). Little is known at present about 
the excitation mechanism. To demonstrate the effect, how- 
ever, we artificially increase the oscillation amplitude ten- 
fold (with respect to the zeus-mp output) by subtracting 
the time average from the signal, multiplying the remaining 
signal by ten, and flnally adding the time average again. The 
results appear in Fig. |5] The sidebands are clearly visible, 
separated by A/ — 17 Hz from the main peak. The width 
of the sidebands is set by the damping rate. 

In principle, a high resolution spectrum of the grav- 
itational wave signal allows us to measure Mg and the 
surface magnetic field fro m A/ jPavne fc Melato^ l2006al : 
IVigelius fc Melatosll2008"al ). A detailed analysis of the global 
MHD oscillations, including the linear response to a stochas- 
tic excitation, will be attempted in a forthcoming paper. 

4.2 Triaxial, nonvibrating mountain 

The Parker instability experienced by an initially axisym- 
metric mountain saturates in a slightly nonaxisymmetric 
state, with ~ Bp and Q12/Q33 ~ 10~^, as discussed in 



10 M. Vigelius et al. 




Figure 8. Fourier transform of tlie gravitational wave signal from a standard neutron star (Ma = 1.4Mc) with an axisymmetric magnetic 
mountain at a distance D = 10 kpc, with wobble angle 9 = 7r/3 (top panels) and d = 0.1 (bottom panels), and inclination i = 7r/3. The 
left (right) panels show the discrete Fourier transforms of the signals /i+(<) [/ix(i)] in the plus (cross) polarisations. The peak is not a i5 
function due to the finite resolution of the discrete Fourier transform. The dashed curves display the detection threshold htjj for Initial 
(upper) and Advanced (lower) LIGO, assuming an observation time of To = 14, a false alarm rate of 1 per cent, and a false dismissal 
rate of 10 per cent. 



see sections 12. II and 12. 21 A third, unequal principal moment 
of inertia i ntroduces new f eatur es into the gravitational wave 
spectrum. IZimmermannI l|l980l ) computed the waveform in 
a small- wobble-angle expansion. IVan Den Broeclj (|2005l ) ex- 
tended this analysis up to second order in 9, finding 



hx {t) = ^ [^x,fc sva{n\kt) + Ax,i, sin(r22fei)] 



(16) 



with ^2*; = 2fiiot + 2fcfiproc and n2fe = ^rot + (1 + 2k)Q,prcc, 

where the fundamental angular frequencies are 



h+{t) = ^ cos(n2fci) + ^+,fc cos(r22fei)] 



(15) 



(h-h){h-Ii) 
hh 



1/2 



(17) 



and 



and 
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Figure 10. Fourier transform of the gravitational wave signal from a standard neutron star {Ma = lAMc) with an nonaxisymmetric 
magnetic mountain at a distance D = 10 kpc, with wobble angle 6 = 0.1 and inclination i = tt/3. The left (right) panels show the 
discrete Fourier transform of the signals h^{t) [hx{t)] in the plus and cross polarisations respectively (second order in 9). The dashed 
curves display h^^ of LIGO and advanced LIGO, assuming an observation time of Tq = 14 d. The lower panels zooms in on the peak at 
2(f2rot + f2prcc), revealing the split peak due to triaxiality. 



^rot — ^ 



1 + - 



i94(i7ra) 



(18) 



In p7|) and H18|) . I\, I2, and /a denote the principal moments 
of inertia, J = {ll^l + /|ni + ll^lY'^ is the (conserved) 
total angular momentum, and f2 = (f2i, SI2, fis) is the angu- 
lar velocity vector. At t = 0, we have fl = (a, 0, b) without 
loss of generality, suc h that a/b ~ 9 for small wobble angles. 
For consistency with IVan Den Broec^ (|2005l ). the principal 
moments of inertia are ordered such that Ii ^ I2 ^ I3 and 
^ 2EI2, where 2E = hnj + I2O.I + I3O.I is the (con- 
served) total kinetic energy. The parameter m is defined by 



{l2-Ii)Iia^ 
{h ~ l2)hb^ ' 



(19) 



K{m) is the complete elliptic integral of the first kind, 
^4 is the fourth Jacobi theta function, with nome q = 
exp[— 7riv(l — m) /K{m)], and q is a solution of 



sn[2iaK (m) , m] 



ihb 
Iia ' 



(20) 



where sn is the (doubly periodic) Jacobi elliptic function. 

The amplitudes A+ and Ax for the plus an d cross po- 
larisations take the form (|Van Den Broec^l2005h 
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Figure 9. Fourier transform of h-|_(t), for a neutron star {Ma = 
lAMc) at a distance D = 10 kpc, with wobble angle d = 0.1 and 
inclination i = 7r/3 (cf Fig.[8ll. The amplitude of the oscillations 
in € has been increased artificially by a factor 10 in order to bring 
out the sidebands. The nonzero width of the sidebands stems from 
numerical damping. 



^fe"(l+COs"i)(/2-/l), 



D 



A\^^ = ^b'[h - (/i + h)/2]-i'\l + cos' i), 



Al a = 



ix,0 

I 

X ,1 
II 

+ ,0 

II 

+,1 

^ V ft — 



- — b'ih - Ii) cosi, 



4 , 
0, 



■-x ,0 
and finally 



(Ji +/2)/2]7'cosi, 
(/i + /2)/2]7sin(2i), 

(/i + J2)/2]7sini, 



(21) 

(22) 

(23) 

(24) 

(25) 
(26) 
(27) 

(28) 

where we omit the factor G/c^ to avoid confusion. In (|21|| - 
(|28p . i denotes the inclination angle, D is the distance to the 
star, and 7 <C 1 is the expansion parameter 

Equations (|15|l and (|16p point to the existence of three 
distinct lines in the gravitational wave spectrum. Line I, 
at 2^10*, stems from the departu r e from axisymmetry [cf. 
the discussion in lVan Den Broeckl (|2005h ]. A body rotating 
around one of its principal axes of inertia looks identical af- 
ter half a period. Line II, at Orot + fiprcc, results from the 
free precession of a nearly axisymmetri c (/i — > /2) object i n 
the small- wobble-angle approximation (|Zimmermannlll980l ) . 
Line III, appearing as a sidelobe to line I at fig + 2r2pi-cc, 



results from the second order expansion of the precession. 
For a biaxial star, lines I and III coincide at the frequency 
Qiot + f2prcc ~ 'in ft,, consistent with equation ((S}. Nonax- 
isymmetry then separates line I, which is shifted to a lower 
frequency and the second-order line III, which remains at 

We compute the gravitational wave spectrum for the 
triaxial equilibrium (Fig.lTJ)) from Qij at t = IOta. In order 



to rescale Qij to a realistic curvature (se e section 12.511, we 
empl oy the relation e = Q33/2I3 oc (jMelatos fc Pavn3 
l2005h . which im plies Q33 oc Mi,Rl, and hen ce Qij oc MtRt 
consistent with IVigelius fc MelatosI (|2008al ) ■ The (up- 
scaled) principal moments of inertia evaluate to [/i/Tr(/)] — 



and 

^2 



1 = -3.9 X 10-^ [/2/Tr(7)] - 1 = -1.8 x 10" 
[J3/Tr(J)]-l = 5.8xl0"^ where Tr(/) = 3.3 xlO'"^ gem' is 
the trace of the moment of inertia tensor. We also consider 
a star with D = 10 kpc, n/2TT = 600 Hz, and 6 ^a/h = 0.1 
rad. This particular choice results in a = 0.93, Slprcc = 0.032 
rad s~^, and f2rot ~ 3.7 x 10^ rad s^^. 

Fig. [To] depicts the spectrum (solid curves) emitted by 
the above star, along with /ith for Initial LIGO and Ad- 
vanced LIGO (dashed curves) , assuming an observation time 
of To = 14 d. Again, all the peaks in this mock spectrum 
are 5 functions in reality; their width stems from the finite 
resolution of the discrete Fourier transform. The peaks at 
fig and Q}q are clearly visible in both polarisations. Neither 
surpasses the noise floor unless the tunability of Advanced 
LIGO is exploited. A larger wobble angle improves matters. 
The size of 9 is controlled be extraneous factors, such as the 
angle CI makes with the accretion disk (Lai 2007), t he accre- 
tion history, the crystallisatio n history of the crust (|Melatosl 



I2OO0I : iMelatos fc Pavn3 



2005^■ and dissipative pr ocesses in 



the superfluid interior (|Cutler fc LindblomI 119871 ). none of 
which are well understood. However, if deeply modulated, 
persistent X-ray pulsations emanate from a hot spot on the 
stellar surface at the mag netic poles, then 6 must be appre- 
ciable (|Chung et al.ll2008l ). 

The bottom panels of Fig. [10] zoom in on the peaks at 
f^o'i revealing the second-order peak at ■ We measure the 
amplitudes to be in the ratios |^+^o/^+,il ~ l^x,o/^x,il ~ 
25. This shows clearly how triaxiality shifts the main peak to 
a lower frequency (2nrot), while only a small, second-order 
peak remains at 2{Q,rot + fiprec), as for a biaxial star. The 
peaks are separated by An/27r = Q.^rac/n = 0.01 Hz. 

Can the extra, "triaxial" peak be exploited to facilitate 
detection, i.e. to increase the signal-to-noise ratio d? The 
generalisation of the results of section IT^ is straightforward, 
as long as we assume that the detector can resolve all three 
spectral lines. After averaging over sky position, inclination, 
and polarisation, we find 



{cF) ^ {dl) + (dh) + {d!u) 

with 



(30) 



By applying the same scaling to all components Qij, we effec- 
tively assume that the three-dimensional structure of the moun- 
tain scales homologously with We find empirically that this 
assumption holds for axisymmetric llPavne &: Melato£l2004l.l2007t) 
and n onaxisymmetric configurations [see IVigelius fc MelatosI 
bOOSal) . and section 4.6]. 
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{h-Iif, 

[/3-(/l+/2)/2]V, 



255h(/iii)D2 



[/3-(/i + /2)/2]V 



(31) 
(32) 
(33) 



with /i,ii = f7o'"/(27r) and /m = fi^'. 

If the canonical star above {Ma = 10"^Mq, D = 10 
kpc, and S — 0.1) is observed with Initial LIGO, we expect 
(df ) = 1.6, (dfi) = 1.1, and (d?„) = 0.003. Importantly, if 
one does not account for the frequency shift of line II when 
searching the LIGO data, one picks up only the strongly 
attenuated line III and hence only half of the total signal- 
to-noise ratio. 

Once an initial detection is made, it is possible to ex- 
tract the inclination, precession angle, deviation from ax- 
isymmetry, and oblate ness parameter 1 — {Ii + I-i)l2Iz 
llVan Den Broeckl [20051 ') from the gravitational wave spec- 
trum. Firstly, i can be found easily from the amplitude ra- 
tio of the two polarisations in any spectral line. Next, the 
amplitude ratio of line II and line III can (knowing i) be 
used to infer Iia/Is,b. The ratio of lines I and II is propor- 
tional to hh{l2 - h)/Iia[h - {h - /i)/2]. The line frequen- 
cies tell us Q,Yot and f^proc- Using the approximate relation 



7r[/3 — (7i -I- I-z) /'2\/2K(m)I-j,, one can finally 



derive the oblateness parameter [/g — (/i + I'2)/2]/Iz- 



4.3 Precession amplitude 

At present, the wobble angle 6 is partly constrained by the- 
ory and observation. We know that it cannot be zero exactly, 
because the fluctuating magnetospheric accretion torque has 
nonzero components perpendicular to f2 (Lai 1999. ), which 
give 6^0 for finite damping. But how small is the steady- 
state 61 If the magnetic axis is aligned with il, due to 
viscous dissipation, before the crust crystallises, then /i, fl, 
and 63 are all aligned. In this case, we see neither preces- 
sion nor X-ray pulsations before accretion begins (provided 
the pulsations stem from a hot spot at the magnetic pole). 
However, a subsequently accreted mountain with a small 
nonaxisymmetry emits gravitational radiation at f — 2f2rot 
(line I), but still no X-ray radiation. On the other hand, if 
the crust crystallises before fl aligns with fi, we see preces- 
sion, pulsation , and gravitational radiat ion in all three lines 
(|Melatosll2000l : iPavne fc Melatosll2006ah . 

Gravitational-wave back-reaction damp s the wobble on 
atime-scale r ~ 10Vr(10"''/e)^(kHz//*)* (|Cutler fc Jones! 
I2OOII ). provided the star can be treated as a fluid body 
with an elastic crust and the precession is torque- free. On 
the other hand, a neutron star with an accreted mountain 
forms a prolate spheroid. In this case, internal dissipation 
in creases the wobble an gl e (9 — > n/2) [cf. the discussion 
in ICutler fc Joned (120011 )]. lAlpar fc SaulisI (|l988l ) examined 
the coupling of the superfluid interior and the stellar crust, 
finding that 6 approaches 7r/2 over ~ 10* free precession 
periods. 

It has been postulated that precession explains the slow 
variation in pulse-arrival-time residuals observed in accret- 
ing neutron stars (e.g. Her X-1) and isolated pulsars (e.g. 
PSR B0959-54 and PSR B1828-11) [cf. the discussion in 



lAkgiin et al] (|2006l )]. IChung et al.] (|2008l ) modelled the X- 
ray fiux modulations observed in XTE J1814-338. They 
found e ^ 10"^ and 0° ^ ^ 5° if t he star precesses 
and e cos 61 ^ 10"^° if not. lAkgiin et al] (2006 1 modelled 
the radio pulses from PSR B1828— 11, finding that they can 
be matched by a precessing, triaxial star, whose degree of 
nonaxisymmetry is small and whose shape is p rolate. The 
wobble angle of t his isolated pulsar is 61 ~ 3° (|Linkll2003l : 
IStah^s et allbOOOl ). 



5 DISCUSSION 

In this article, we give improved estimates for the strength 
and spectrum of gravitational waves from accreting millisec- 
ond pulsars. For the first time, we include the effects of non- 
axisymmetry and resistive relaxation on the gravitational 
wave signal. Furthermore, we justi fy the curvature down- 
scaling introduced in previous w ork (|Pavne fc Melatosll200j . 
l2007l : IVigelius fc Melatoll2008al ) and achieve self-consistent 
mountain configurations with Ma ~ 1O~^M0, ten times 
larger than constructed previously, by growing them from 
scratch by injection. We cannot find any evidence for a grow- 
ing, axisymmetric instability even at these Ma- 

Taken at face value. Figs. |6] and [7] imply that moun- 
tains with (/./O.3kHz)2(M„/lO"'^M0) > 1 are detectable 



in principle even with Initial LIGO. With Advanced LIGO, 
chances are even better over the whole 0.1 — 1 kHz band for 
Ma > 1O ~^M0. While this e ncour aging result was foreshad- 
owed by iMelatos fc Pavnel (|2005l ). this paper includes for 
the first time the dynamics of nonaxisymmetry and resis- 
tive relaxation, which act to weaken the signal - yet still 
the prospects remain bright. In principle, the noise floor 
of Advanced LIGO can be lowered even further by using 
a narrowband configurati on exploiting a squeezed vacuum 
(jBuonanno fc Chenl l200j ). improving the sensitivity by as 
much as threefold for / < 400 Hz. A successful detection 
will let us test the reciprocal dependence of the gravitational 
wave amplitude on the ma gnetic dipole moment p redicted 
for magnetic mountains by Melatos fc Pavng (|2005l ). 

Of course, magnetic mountains have not been detected 
during recen t searches for low-mass X-ray binaries i n the S5 
LIGO data (| Abbott et al.ll2007l : IWatts et al.ll2008l '). More- 
over, for many LMXBs, the quadrupole moment predicted 
by Fig. [5] is too large to be consistent with the mea- 
sured spin frequency, because the accretion-driven recycling 
process is "stalled" by the gravitational-wave spin-down 
torque (|Bildstenlll998l : IChakrabartv et al.ll2003t ). This sug- 
gests that other physical processes exist, unaccounted for 
so far, which reduce the mountain. The circle of candidates 
has shru nk significantly since the ex clusion of resistive re- 
laxatiorFI (|Vigelius fc MelatosI l2008bh , but there are others, 
principally hydrodynamic sinking. In all the modelling to 
date, we assume that the mountain sits on a hard surface. 
Realistically, the neutron star crust i s not impenetrable, and 
part of the mountain sinks into it (|Konar fc Bhattacharval 



^ It is possible that the electrical resistivity is substantially 
higher than the cont ribution from electron-phonon scattering 
llCumming et al.ll2004l) or that the H all effect plays a dominant 
relaxing role llPons fc GeppertJ 1200711 . These possibilities will be 
pursued in future investigations. 
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I1997I : IChoudhuri fc Konaij |2002| : IPavne fc MelatosI l2004l fl. 
An attempt to treat sinking by growing mountains on top 
of a soft crust, modelled by a polytropic equation of state, 
will be presented elsewhere [Wette et al. (in preparation)]. A 
more sophisticated model would include the realistic, strati- 
fied nuclear composition of the cru st and its evolutiori in re- 
sponse to pycnonuclear reactions (|Brown fc Bildstenlll998l : 
lUshomirskv et al.ll20od ). The Coriolis force may also push 
the mountain to wander across the surface, especially in 
the fastest spinning LMXBs, although the enhanced mag- 
netic field in the equat orial belt opposes the wandering 
(jPavne fc MelatosI [2006b[ ). 

For simplicity, we assume an isothermal equation of 
state throughout this article. During the late stages of ac- 
cretion (Ma > lO'^M©), however, the magnetic mountain 
mass is comparable to the mass of the neutron star crust and 
the equilibrium exhibits a wide range of density and temper- 
ature. Pycnonuclear reactions in the deep regions (p > lO'^^ 
g cm~^) feed thermal energy into an adiabatic mountain. 
The assumption of isothermality breaks down and a real- 
istic equation of state for non-catalyzed matter is required 
iHaensel fc Zdunik|[l990ah . In particular, the accreted ma- 
terial is expected to soli dify at densities > 10* g cm""' 
(|Haensel fc Zdunik|[l99qbl) and the crust needs to be mod- 
elled as an elastic solid " fUshomirskv et al.llioool '). The effect 
of a realistic equation of state is subject of current work and 
the results will be presented elsewhere. 

Finally, we note that the full (discrete and continu- 
ous) spectrum of global MHD mountain oscillations contains 
valuable information about the structure of the star, e.g. 
the strength of the surface magnetic field. In principle, this 
spectral information will be accessible by third-generation 
gravitational- wave interferometers with improved sensitivity 
and frequency resolution. As a first step, we will compute 
the continuous part of the MHD spectrum in a forthcoming 
paper. 
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APPENDIX A: GROWING MOUNTAINS WITH 
ZEUS-MP 

In this appendix, we describe one way to grow a moun- 
tain with Mg » M c, using the ideal-MHD code zeus-mp 
I Haves et all 120061) extend ed to include Ohmic dissipation 
I Vigelius fc Melatos|[2008bl l. 

We start with an isothermal atmosphere (initial mass 
<C Mc) resting on a hard surface ai r = R, threaded by 
a dipolar magnetic field. We then inject matter quasistati- 
cally into the atmosphere from below (at r = i?*) along a 
polar subset of the magnetic field lines, imitating disk-fed 
accretion. We artificially boost the mass flux so as to grow 
the mountain in a reasonable time, albeit slowly compared 
to the hydromagnetic Alfven time-scale. This is achieved by 
boosting the injection mass density po while keeping the in- 
jection speed vo below the gravitational escape speed (i.e. 
escape from the simulation volume). 

Mass is added at r = i?, at a rate M{6) = 
povo ex p(— bsin^ ff), where b is the polar cap radius scaled 
to -R, (|Pavne fc Melatod [2003 ). By injecting at r = i?. 
(from "below") instead of r = Rm (from "above"), we ex- 
ploit the flux-freezing property of ideal MHD to avoid having 
to locate the changing intersection point of the flux surface 
'I/ = v&o with r = Rm; cf. the bootstrapping algorithm in 
iPavne fc MelatosI (|2007l ). We find empirically that a good 
compromise between shortening the run time while main- 
taining quasistatic injection is po = 5.1 x 10^" g cm^'^ and 
Vo = 10* cm s"^. 

Although the magnitude of i;o is constant, its direction 
must be parallel to the dipolar magnetic field everywhere at 
r = i?, , in order to respect the flux-freezing constraint. The 
6 — n/2 surface is reflecting [ojb.no js(l)= 5], with nor- 
mal magnetic field, which translates to vx = By =0. The 
line 6 = is also reflecting [ijb.nijs(l)= -1], with tangen- 
tial magnetic field (vx ~ Bx =0). The oute r r surf ace is 
outflow [cf. section 3.2 in I Vigelius fc MelatosI (|2008al )]. 

We first verify that the growing mountain upholds the 
correct mass/flux ratio. Fig. lAll displays dM/dtp for the 
injection simulation (squares), the Grad-Shafranov equilib- 
rium (diamon ds), and the analytic d istribution dM/d%p = 
exp(--i/'/V'a) (|Pavne fc MelatosI |2003 ). All match to better 
than one per cent. While the snapshot is taken at f = 2ta 
(which translates to Ma = 0.46 x lO'^M©), we verify that 
dM/dtf) remains unchanged with time. 

By way of verification, we compare a grown mountain 
to a Grad-Shafranov equilibrium for Ma = 1.2 x 10"* in 
Fig. IA2I The left panels display the same density contours 
for the grown mountain (dashed curves) and the equilibrium 
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Figure Al. dM/dip in arbitrary units as a function of ^/tpa 
for the Grad-Shafranov equilibrium (diamonds) and the injection 
simulation (square). For comparison, we also include the theoret- 
ical value dM/dij) = e^p(—ip/'il)a) (solid curve). The snapshot is 
taken at i = 2tj^, although the mass-flux ratio remains constant 
with time, of course. 
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Figure A2. Meridional section of density contours (left) and 
magnetic field lines (right) for a grown mountain (dashed curve) 
and a Grad-Shafranov equilibrium (solid curve) with Ma = O.lMc 
(upper panels) and Ma = Mc (lower panels) . While the agreement 
in the density is good, discrepancies in the magnetic field arise 
due to the different boundary conditions at r = Rm (see text). 

(solid curves), for Ma = O.lMc (top panels) and Ma ~ Mc 
(bottom panels). The density contours in the region close 
to ^ X — (r — _R*)//io ^ 2 match reasonably well, with 
a deviation of less than one per cent radially and five per 
cent laterally. While the magnetic field lines match close to 
the stellar surface, there is a considerable discrepancy in the 
(less important) outer regions, where B is weak (and indeed 
p is small). This stems from the different implementation 
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of the outer boundary condition at r = Rm- The outflow 
boundary condition in zeus- mp enforces a van i shing gra- 
dient of B at r = _R™, while IPavne fc Melatod lj2004h en- 
force dil)/dr = 0, thereby imposing the additional constraint 
Bg = 0. The exact form of the outer boundary condition 
depends on the interaction of the magnetosphere with the 
accretion disk and is poorly known. We are not concerned 
with the details and just note that the magnetic field decays 
oc r~^, such that it is several orders of magnitude lower than 
at the surface. 



